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Abstract 

This  work  compares  the  performance  of  lithium  batteries  with  polymer  electrolytes  with  unity  (“ionomer”)  and  nonunity  (‘"polymer 
electrolyte”)  transference  numbers.  The  study  is  performed  with  respect  to  a  particular  cell  chemistry,  Li  metal  |  polymer  |  Li V6013- 
composite  electrode,  which  is  currently  a  top  candidate  for  use  in  electric  vehicles.  Cell  performance  was  modeled  to  determine  the  best 
possible  performance  of  cells  containing  four  different  electrolytes:  “ideal”  polymer  membrane  and  ionomer  with  properties  defined  by 
USABC  goals,  and  the  presently  best  available  polymer  electrolyte  and  ionomer.  Positive  electrode  thickness,  porosity,  and  current 
density  were  varied  to  find  the  cell  geometry  with  the  highest  combined  energy  density  and  peak  power  performance  for  cells  with  each 
electrolyte,  and  concentration  and  potential  profiles  are  examined  to  determine  the  limitations  of  the  electrolytes.  The  results  show  that  at 
40°C,  the  “ideal”  polymer  electrolyte  can  provide  104  W  h/kg  and  99  Wp/kg,  the  “ideal”  ionomer  can  provide  94  W  h/kg  and  58 
Wp/kg,  and  the  currently  available  electrolytes  can  provide  about  one-fifth  of  these  values.  Published  by  Elsevier  Science  S.A. 
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1.  Introduction 

This  study  analyzes  electrochemical  performance  of 
polymer  electrolytes  in  lithium-polymer  batteries,  the 
trade-offs  between  unity  (“ionomer”)  and  nonunity 
(“polymer  electrolyte”)  transference  numbers,  and  exam¬ 
ines  the  basis  for  long-term  performance  goals  for  polymer 
electrolytes  for  use  in  electric  vehicles  that  have  been 
developed  by  the  United  States  Advanced  Battery  Consor¬ 
tium  (USABC).  The  analysis  is  performed  with  respect  to 
a  particular  cell  chemistry,  Li  metal  |  polymer  |  LiV6013- 
composite  electrode,  which  is  currently  a  top  candidate  for 
use  in  electric  vehicles.  Cell  performance  was  modeled  to 
determine  the  best  possible  performance  of  cells  contain- 
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ing  four  different  electrolytes:  “ideal”  polymer  electrolyte 
and  an  “ideal”  ionomer  with  properties  defined  by 
USABC  goals,  and  the  presently  “best-available”  (best 
of  those  tested  in  this  laboratory)  polymer  electrolyte, 
oxymethylene-linked  poly(ethylene  glycol)  (PEMO) 
(molecular  weight  400)  with  lithium  bis(trifluoromethyl- 
sulfonyl)  imide  (TFSI),  and  a  best  available  ionomer,  a 
random  copolymer  of  methylpolyethylene  glycol  acrylate 
(molecular  weight  500)  and  lithium  sulfonated  diethylene 
glycol  acrylate  (see  Table  1).  The  transport  properties  of 
an  ionomer  are  completely  described  by  the  ionic  conduc¬ 
tivity,  which  was  measured  using  the  technique  described 
in  Ref.  [5].  Transport  in  a  polymer  electrolyte  involves  two 
additional  properties.  Appendix  A  describes  the  complete 
transport  properties  of  the  polymer  electrolytes  used  in  the 
simulations. 

The  positive  electrode  thickness,  porosity,  and  current 
density  were  varied  to  find  the  combination  which  yields 
the  highest  product  of  specific  energy  for  a  3-h  galvanos- 
tatic  discharge  and  peak  specific  power  for  30  s  after  80% 
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Table  1 

Transport  properties  (ionic  conductivity,  cation  transference  number,  and 
salt  diffusion  coefficient)  of  the  electrolytes  modeled  in  this  study 

Ideal  case  Best  currently  available 

(USABC  goals) 


Polymer  electrolyte 


k  (S/cm) 

10“3 

10~4 

t° 

+ 

0.3 

0.1 

D  (m2/s) 

10“n 

6X10 

Ionomer 

k  (S/cm) 

10“4 

4X10 

of  the  discharge  time  (at  2.4  h).1  The  model  used  has  been 
described  by  Fuller  et  al.  [1].  The  mass  used  to  calculate 
specific  energy  and  peak  specific  power  includes  the  mass 
of  the  current  collectors,  active  material,  conductive  filler, 
and  electrolyte,  but  does  not  include  packaging.  The  model 
was  run  for  values  of  current  density  ranging  from  0.05  to 
1.7  mA/cm2,  positive  electrode  thickness  ranging  from  10 
to  160  ptm,  and  electrolyte  volume  fraction  in  the  positive 
electrode  ranging  from  0.15  to  0.5.  For  a  given  application, 
the  required  total  current  will  be  determined  by  the  load. 
The  current  density  required  for  this  load,  then,  is  the  total 
current  divided  by  the  cell  area.  Thus,  optimizing  the 
current  density  allows  one  to  calculate  the  cell  area  re¬ 
quired  for  a  given  load.  All  other  parameters  used  in  the 
model  are  shown  in  Appendix  B. 

2.  Results  and  discussion 

Table  2  shows  the  results  of  the  simulations  for  an 
operating  temperature  of  40°C.  One  should  note  that  the 
maximum  peak  specific  power  and  specific  energy  listed 
are  those  that  give  the  highest  combination  of  energy  and 
power:  cells  can  be  designed  to  increase  one  at  the  expense 
of  the  other,  within  the  capabilities  of  the  cell  chemistry, 
as  needed  by  a  particular  application.  The  “percent  derat¬ 
ing”  can  be  defined  as  the  ratio  of  the  maximum  specific 
energy  when  the  cell  is  optimized  for  specific  energy  to 
that  when  the  cell  is  optimized  for  specific  energy  and 
peak  specific  power  combined.  For  these  cells,  the  derating 
is  approximately  50%. 

The  polymer  electrolyte  cells  can  provide  greater  spe¬ 
cific  power  than  the  ionomers  because  of  their  higher  ionic 
conductivity.  Performance  of  the  ionomers  is  limited  by 
the  ohmic  resistance,  while  the  polymer  electrolytes  are 


1  Optimizing  the  product  of  specific  energy  and  peak  specific  power 
yields  a  result  very  close  to  that  obtained  when  the  cell  is  optimized  by 
finding  the  largest  point  which  lies  on  a  line  whose  slope  is  the  ratio  of 
the  maximum  specific  energy  to  maximum  peak  specific  power,  when  the 
cell  is  optimized  for  each  performance  criterion  separately,  as  described 
in  Ref.  [2]. 


Table  2 


Best  combined  energy-power  performance  of  Li  |  LiV6013  cells  with  four 
electrolyte  materials,  with  50  |xm- thick  separators 


Polymer  electrolytes  Ionomers 

Available 

Ideal 

Available 

Ideal 

Specific  energy  (W  h /kg) 

24 

104 

23 

94 

Peak  specific  power  (W /kg) 

20 

88 

8 

58 

Current  density  (mA/cm2) 

0.15 

1.0 

0.15 

0.95 

Positive  electrode  thickness 
(|xm) 

20 

90 

20 

90 

Porosity 

0.3 

0.4 

0.4 

0.35 

limited  by  the  formation  of  concentration  gradients  that 
lead  to  salt  depletion  at  the  cathode  and  risk  of  precipita¬ 
tion  at  the  anode.  The  ideal  electrolytes  can  be  used  in 
thicker  cells  with  more  active  material  per  unit  area  be¬ 
cause  transport  to  the  active  material  is  not  hindered  as 
much  as  in  the  best-available  electrolytes. 

One  can  gain  much  insight  into  the  cell  performance  by 
examining  the  state  of  the  cell  at  the  end  of  the  3-h 
discharge.  For  the  ionomer,  the  electrolyte  concentration  is 
constant.  This  restriction  has  two  major  effects.  First,  the 
reaction-rate  distribution  in  the  positive  electrode  is  nearly 
uniform,  varying  only  due  to  ohmic  drop  in  the  electrolyte. 
Therefore,  the  active  material  (Fig.  1)  is  filled  nearly 
uniformly  across  the  electrode.  Second,  there  is  no  concen¬ 
tration  overpotential.  As  a  result,  the  cell  with  the  highest 
specific  energy  will  be  one  that  has  just  enough  active 
material  to  supply  the  load  for  the  design  time  of  3  h,  after 
which  the  cell  is  completely  discharged  and  the  cell  volt¬ 
age  will  fall  abruptly.  The  currently  available  ionomer  has 
high  ohmic  resistance  (Fig.  2),  and  therefore  less  of  the 


Fig.  1.  Active  material  utilization  (expressed  as  y  in  LivV6013,  where  y 
can  range  from  0  to  8)  across  the  positive  electrode  for  the  best-available 
and  ideal  ionomers,  with  50  |xm  separator  thickness,  after  2.4  and  3  h  of 
discharge  at  the  optimum  current  density  listed  in  Table  2.  The  active 
material  reacts  nearly  uniformly  due  to  the  constant  electrolyte  concentra¬ 
tion.  A  higher  fraction  of  the  active  material  can  be  utilized  with  the  ideal 
electrolyte  before  the  cutoff  potential  is  reached  because  of  the  lower 
ohmic  losses. 
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Fig.  2.  Electrolyte-phase  potential  in  the  best-available  and  ideal  ionomers,  after  2.4  and  3  h  of  discharge  at  the  optimum  current  density  listed  in  Table  2, 
and  potential  at  the  current  collectors,  demonstrating  the  ohmic  resistance  in  the  electrolytes. 


active  material  is  utilized  before  the  cutoff  voltage  is 
reached.  Thinner  electrodes  minimize  this  ohmic  loss. 

The  performance  of  the  polymer  electrolytes  is  domi¬ 
nated  by  the  concentration  profiles  (Fig.  3).  If  the  concen¬ 
tration  in  the  cathode  is  driven  to  zero,  the  cell  voltage 
quickly  drops.  The  cell  dimensions  and  current  density  for 
any  given  application  must  be  designed  to  maintain  the  salt 
concentration  in  the  cathode  and  to  prevent  salt  precipita¬ 
tion  at  the  anode.  Concentration  gradients  in  the  PEMO- 
TFSI  are  exacerbated  by  the  low  transference  number  (see 
Fig.  6).  Therefore,  only  low  current  densities  can  be  used 
in  these  cells,  and  most  of  the  active  material  is  not 
utilized  (Fig.  4).  In  contrast,  over  75%  of  the  capacity  of 
the  cell  with  the  ideal  polymer  electrolyte  can  be  used.  Fig. 
5  shows  that  the  PEMO-TFSI  cell  has  a  potential  above 
the  cutoff  potential  at  the  end  of  discharge.  If  one  were  to 
design  a  cell  to  maximize  energy  density  for  a  3-h  dis¬ 
charge,  the  cell  would  just  reach  the  cutoff  potential  at  the 
end  of  3  h.  However,  in  this  study,  we  have  optimized  the 
cells  for  the  highest  product  of  specific  energy  and  specific 


Fig.  3.  Electrolyte  concentration  profile  across  the  separator  and  positive 
electrode  in  the  best-available  and  ideal  polymer  electrolytes  after  2.4  and 
3  h  of  discharge  at  the  optimum  current  density  listed  in  Table  2. 


power.  Higher  current  densities  cause  larger  concentration 
gradients  that  greatly  inhibit  the  cell’s  power  pulse  perfor¬ 
mance.  In  order  to  maintain  a  reserve  for  a  power  pulse, 
lower  current  densities  must  be  used. 

USABC  goals  specify  a  current  density  of  1  mA/cm2 
for  a  design  discharge  time  of  3  h.  This  current  density 
does  provide  the  maximum  combination  of  specific  energy 
and  peak  specific  power  for  the  ideal  electrolytes.  Higher 
current  densities  would  provide  higher  specific  energy  at 
the  cost  of  peak  power  performance,  and  conversely  lower 
current  densities  would  give  higher  peak  specific  power 
and  lower  utilization  of  the  active  material.  The  best 
currently  available  electrolyte  materials  cannot  supply  such 
high  current  densities.  The  ohmic  drop  across  the  low-con¬ 
ductivity  ion  exchange  membrane  would  drop  the  voltage 


Fig.  4.  Active  material  utilization  across  the  positive  electrode  for  the 
best-available  and  ideal  polymer  electrolytes  after  2.4  and  3  h  of  dis¬ 
charge  at  the  optimum  current  density  listed  in  Table  2.  Because  of  the 
large  concentration  gradients  shown  in  Fig.  3,  the  reaction  rate  is  highly 
nonuniform  and  much  of  the  electrode  is  not  utilized. 
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Fig.  5.  Electrolyte-phase  potential  in  the  50  |xm  best-available  and  ideal 
polymer  electrolytes,  after  2.4  and  3  h  of  discharge  at  the  optimum 
current  density  listed  in  Table  2,  and  potential  at  the  current  collectors. 


below  tolerable  limits  for  electronic  needs.  Concentration 
gradients  in  the  polymer  electrolyte  would  cause  a  risk  of 
salt  precipitation  and  a  high  concentration  overpotential. 

The  conductivity  at  40°C  of  PEMO-TFSI  at  its  maxi¬ 
mum  is  the  same  as  that  of  the  ideal  ionomer.  However, 
the  performance  is  clearly  worse,  due  to  the  effects  of 
concentration  gradients.  In  addition,  while  the  conductivity 
of  the  ideal  ionomer  is  an  order  of  magnitude  lower  than 
that  of  the  ideal  polymer  electrolyte,  its  specific  energy 
performance  is  only  10%  lower.  From  these  observations, 
we  conclude  that  the  most  significant  improvements  in 
polymer  electrolytes  would  be  obtained  from  an  increase 
in  the  transference  number. 

Remarkably,  at  85°C,  the  properties  of  the  electrolytes 
investigated  in  this  work  are  very  close  to  the  USABC 
goals:  /+  is  slightly  higher  (0.5  instead  of  0.3),  D  slightly 
lower  (5  X  10“ 12  instead  of  1  X  10~u  nr/s,  and  the 
conductivity  at  its  maximum  is  the  same  value  (1  X  10“ 3 
S/cm).  Therefore,  the  performance  of  the  polymer  elec¬ 
trolyte  at  85 °C  will  be  similar  to  the  performance  of  the 
ideal  polymer  electrolyte.  The  question  then  arises  as  to 

Table  3 

Best  combined  energy-power  performance  of  Li  |  LiV6013  cells  with 
available  polymer  electrolyte  and  ionomer  electrolyte  materials,  with  10 


|xm-thick  separators 

Available  Available 

polymer  ionomer 

electrolyte 

Specific  energy  (W  h/kg)  37  41 

Peak  specific  power  (W /kg)  23  20 

Current  density  (mA/cm2)  0.20  0.25 

Positive  electrode  thickness  (|xm)  20  30 

Porosity  0.4  0.4 


Fig.  6.  Transference  number  of  LiTFSI  in  PEMO  at  40°C. 


whether  it  is  better  to  focus  on  designing  an  application 
and  cell  that  can  operate  with  presently  available  polymer 
electrolytes  and  at  85°C,  or  to  focus  on  synthesizing  better 
polymer  electrolytes  that  can  meet  the  USABC  goals  at 
room  temperature. 

For  this  study,  50  pirn  was  selected  as  the  separator 
thickness,  in  order  to  maintain  a  safety  margin  to  ensure 
separation  of  the  electrodes.  Future  manufacturing  methods 
may  allow  thinner  separators  to  be  used.  To  investigate  the 
performance  enhancements  possible  from  simply  making 
the  separator  thinner,  cells  were  optimized  for  the  best 
currently  available  polymer  electrolyte  and  ionomer  with  a 
separator  thickness  of  10  pum.  Results  are  shown  in  Table 
3.  For  separator  of  one-fifth  the  thickness,  we  see  that  the 
specific  energy  almost  doubles.  The  ionomer  is  able  to 
deliver  over  twice  as  much  power  because  of  the  decreased 
ohmic  drop  across  the  separator  (note,  however,  that  the 
power  does  not  increase  by  five  because  much  of  the 
resistance  is  in  the  porous  positive  electrode).  The  peak- 
power  output  of  the  polymer  electrolyte  does  not  change 
significantly  because  this  configuration  favored  specific 
energy  over  peak  specific  power.  We  see  that  while  thinner 
separators  do  result  in  improved  energy  density,  the  in¬ 
crease  is  not  simply  proportional  to  the  decrease  in  thick¬ 
ness,  and  greater  performance  improvements  would  result 
from  improved  transport  properties. 


Fig.  7.  Salt  diffusion  coefficient  of  LiTFSI  in  PEMO  at  40°C. 
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Fig.  8.  Log(ionic  conductivity)  of  LiTFSI  in  PEMO  at  40°C. 

We  note  that  positive  electrode  chemistries  with  higher 
potentials  may  be  selected  for  use  in  lithium-polymer  cells. 
The  chemistry  selected  for  this  study  is  that  described  by 
West  et  al.  [3]  The  shape  of  the  open-circuit  potential 
profile  (i.e.,  the  shape  of  sloped  and  plateau  regions)  has  a 
significant  effect  on  cell  performance,  particularly  in  the 
polymer  electrolytes,  due  to  the  effect  on  the  reaction  rate 
distribution  across  the  positive  electrode.  Therefore,  a  sep¬ 
arate  study  would  be  necessary  to  assess  accurately  the 
performance  of  chemistries  with  differently  shaped  open- 
circuit  potential  profiles.  The  effect  of  merely  increasing 
the  open-circuit  potential  without  changing  its  shape  on  the 
specific  energy  of  a  cell  with  a  given  geometry  is  linear; 
i.e.,  the  change  in  specific  energy  that  would  result  from 
using  a  positive-electrode  chemistry  with  an  open-circuit 
potential  1  V  higher  than  that  used  in  this  study  would  be 
(1  V)  X  (current  density)  X  (discharge  time)/(cell  mass 
per  unit  area),  or  about  9  W  h/kg  for  the  best  available 
polymer  electrolyte  examined  here.  The  effect  on  peak 
specific  power  is  more  complex  and  must  be  determined 
by  separate  simulations;  for  PEMO-TFSI,  the  increase  in 
peak  specific  power  due  to  an  increase  in  cell  potential  of 


Fig.  9.  Electronic  conductivity  in  S/m  as  function  of  state  of  charge, 
from  Ref.  [3]  corrected  for  porosity. 


1  V  is  approximately  15  W/kg.  Note  that  the  above 
discussion  assumes  a  fixed  cell  geometry  and  fixed  differ¬ 
ence  between  the  initial  potential  and  cutoff  potential. 

3.  Conclusion 

In  this  analysis,  we  have  shown  that  cells  operated  at 
40°C,  using  the  best  polymer  electrolyte  and  ionomer 
synthesized  to  date  in  our  laboratories,  would  have  approx¬ 
imately  one-fifth  the  specific  energy  and  peak  specific 
power  of  cells  using  electrolytes  that  met  US  ABC  goals. 
The  optimum  geometry  of  porous  electrodes  depends 
strongly  on  what  one  is  optimizing  for  (i.e.,  energy  or 
power)  and  the  transport  properties  of  the  electrolyte. 
Rather  than  specifying  a  current  density,  one  should  opti¬ 
mize  for  the  separator  area  that  gives  the  best  overall  cell 
performance.  The  polymer  electrolytes  have  severe  trans¬ 
port  limitations.  Research  into  polymer  electrolytes  should 
focus  on  increasing  the  diffusion  coefficient  and  transfer¬ 
ence  number,  not  just  the  conductivity. 

The  direction  of  future  research  efforts  should  be  guided 
by  the  needs  of  the  intended  applications.  Further  simula¬ 
tions  are  necessary  to  determine  whether  the  available 
electrolytes  could  fulfill  the  energy  demands  of  current 
applications  such  as  laptop  computers  or  electric  vehicles. 
There  is  a  tradeoff  between  designing  the  application  and 
cell  to  operate  at  higher  temperatures  versus  trying  to 
synthesize  polymers  with  better  transport  properties. 
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Appendix  A.  Transport  properties  of  polymer  elec¬ 
trolytes  used  in  the  simulations 

The  ideal  polymer  electrolyte  was  assumed  to  have 
transport  properties  which  are  independent  of  concentra¬ 
tion.  The  values  used  were  an  ionic  conductivity  of  0.1 


S/m  and  transference  number  of  0.3,  from  USABC  goals; 
a  salt  diffusion  coefficient  of  1  X  10_u  m2/s;  and  an 
assumed  thermodynamic  factor  of  1.0. 

Properties  of  PEMO-TFSI  were  measured  at  40°C, 
using  the  methods  outlined  in  Refs.  [4,5].  An  average  value 
of  the  thermodynamic  factor,  which  is  defined  as  1  + 
d(ln/)/d(lnc),  where  /  is  the  mean  ionic  activity  coeffi¬ 
cient,  of  4.3  was  used.  Least-squares  fits  to  the  concentra¬ 
tion  dependence  of  the  transference  number,  salt  diffusion 
coefficient,  and  ionic  conductivity  are  shown  in  Figs.  6,  7 
and  8,  respectively. 


Appendix  B.  Parameter  values  for  simulation  of  Li  metal  I  separator  |  LiV6013 


Parameter 
Li  foil  thickness 

Separator  thickness 

Current  collector  thicknesses 

Temperature 

Initial  salt  concentration 

Initial  y  in  Li  vV6Ol3 
Cutoff  potential 

Solid  diffusion  coefficient  in 

positive  material 

Positive  material  particle  radius 

Volume  fraction  of  carbon  in  cathode 

Electronic  conductivity  of  composite  positive  electrode 

Exchange  current  density  at  negative  electrode 

Exchange  current  density  at  positive  electrode 

Film  resistance  at  negative  electrode 

Coulombic  capacity  of  negative  electrode 

Coulombic  capacity  of  positive  electrode 

Density  of  negative  active  material 

Density  of  positive  active  material 

Density  of  carbon  filler 

Density  of  polymer  electrolyte 

Density  of  copper  current  collector 

Density  of  aluminum  current  collector 

Open-circuit  potential 

Transfer  coefficient 


Value 

Explanation /reference 

50  pm 

Capacity  of  negative  =  2.5  times  capacity 
of  positive  electrode 

Minimize  thickness  while  preventing 

25  pm 

short-circuiting 

Standard  size  (1  mil)  for  Al  and  Cu  foil 

313  K 

Just  above  room  temperature 

1300  mol/m3 

Conductivity  maximum  for 

0.3 

polymer  electrolytes 

Fully  charged  [3] 

1.9  V 

Prevent  lattice  degradation, 

1  X  1(T13  m2/s 

protect  electronics  [3,6] 

[3,6] 

10  pm 

Estimate  to  balance  power  performance 

0.2 

with  degradation  reactions 

See  note  1 

11.45  S/m 

See  note  1 

3.1  mA/cm 2 

[7].  Reference  concentration  is  1  M 

0.4  mA/cm2 

[3].  Reference  concentration  is  1  M  and  y 

0.01  D  m2 

Based  on  Ref.  [8] 

3862.5  mA  h/g 
417.4  mA  h/g 

Based  on  y  in  L/V6013  ranging  from  0  tc 

534  kg/m3 

[9] 

3900  kg/m3 

Interpolated  from  Ref.  [9] 

1800  kg/m3 

[9] 

1700  kg/m3 

Measured 

8954  kg/m3 

[9] 

2707  kg/m3 

[9] 

0.5 

See  note  2 
[7] 

B.l.  Note  1 

West  et  al.  [3]  measured  the  electrical  conductivity  of  a 
porous  compressed  electrode  of  Li}V6012  864  as  a  function 


of  y.  The  electrical  conductivity  ranged  from  40  S /m  at 
y  =  0.6  to  0.4  S /m  at  y  =  6.6,  with  an  average  of  7  S /m 
(corrected  for  a  porosity  of  about  0.6).  Fig.  9  shows  a  fit  to 
the  data  from  Ref.  [3],  using  the  equation:  cr  =  9  + 
120e(-L5;y)  —  9tanh(  y  —  4.5). 
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Using  the  model  of  Meredith  and  Tobias  [11]  for  a 
concentrated  mixture  of  spheres  of  carbon  filler  with  a 
modest  particle  size  distribution  in  a  matrix  of  vanadium 
oxide,  it  is  calculated  that  a  20%  volume  fraction  of  carbon 
yields  a  minimum  effective  conductivity  of  0.64  S /m  and 
an  average  conductivity  of  1 1 .45  S /m.  The  value  used  for 
the  electrical  conductivity  of  carbon  was  105  S/m  [12]. 
The  Bruggeman  correction  is  then  used  to  correct  the 
effective  conductivity  of  the  matrix  for  its  porosity. 

B.2.  Note  2 

The  open-circuit  potential,  U,  of  LiyV6013  and  related 
nonstoichiometric  vanadium  oxides  has  been  reported  by 
Refs.  [3,6,10],  all  of  which  report  very  similar  profiles 
(Fig.  10).  The  equation  used  in  this  simulation  was  ob¬ 
tained  by  a  fit  to  data  from  Ref.  [3]  at  25 °C,  for  y  ranging 
from  0  to  8.25: 

U=  1.9  +  0.13  tanh  (  —  1 .67y  +  1.7) 

+  0.2  tanh  (  —  1.67y  +  6.2) 

+  0.56  tanh  (  —  3.33  y  +  29.5). 
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